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Figure 2: Yi Wan, Gerardo Ortiz and Philip Phillips 




10 ' ' ' ' ' ' ' ' ' ' ' 

2 4 6 8 10 

B [Tesia] 



Figure 3: Yi Wan, Gerardo Ortiz and Philip Phillips 



Pair Tunneling in Semiconductor Quantum Dots 



Yi Wan, Gerardo Ortiz, and Philip Phillips 
Loomis Laboratory of Physics 
University of Illinois at Urbana- Champaign 
1100 W.Green St., Urbana, IL, 61801-3080 

We propose here a model for the pair tunneling states observed by Ashoori and co-workers (Phys. 
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state. The physical potential in which the two electrons are bound can be composed of a Si impurity 
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Single electron capacitance spectroscopy-'^ (SECS) has 
made it possible to detect single electron charging events 
in a semiconductor quantum dot that otherwise con- 
tains no excess electrons. In this regime, the Coulomb- 
blockade energy is negligible relative to the energy of local 
potential fluctuations. As the bias voltage is increased, 
electrons tunnel into local defect states on the quantum 
dot. Beyond a critical value of the bias voltage, sufficient 
charge has accumulated on the dot to form a uniform 
charging state. Transport becomes Coulomb-limited for 
all subsequent charging events^'^. Recently, Ashoori and 
CO- workers"' reported a series of measurements on a GaAs 
tunnel capacitor in the sub-Coulomb blockade regime. 
A tunnel capacitor is a layered structure consisting of a 
source lead, a tunnel junction, a laterally-confined GaAs 
quantum well, a (5— doped backing layer of AlGaAs and 
a gate electrode^ . The (5-dopants in the backing layer 
ionize creating an excess of negative charge that repels 
the tunneling electrons from the backing layer. As a con- 
sequence, charge is restricted to fiow back and forth be- 
tween the source lead and the quantum well. The novel 
feature of these experiments is that single as well as two- 
electron tunneling events were observed^ in quantum dots 
Ijim in diameter. Direct confirmation of pair (or two 
electron)-tunneling is obtained from the evolution of the 
tunneling peaks in the presence of a magnetic field. The 
pair peaks, labeled X, Y , and Z in Fig. 6b of Ref. [lb], 
were shifted to higher gate voltages (Vq) and were ob- 
served to split into two distinct tunneling peaks at critical 
values of the field. Be 2T. Above Bg, the spectroscopic 
signal now corresponds to two energetically-resolved tun- 
neling events. The amplitude for each peak is consistent 
with that for single-electron tunneling. Hence, above Be, 
each state carries a single electron, whereas the charge 
in the states X, Y, and Z at zero-field must be 2e~. At 
a field of 10 T, the splitting between the two states for 
all three pair states is roughly 3 meV. The Zeeman en- 
ergy at this field strength is 0.25 meV. Hence, breakup of 
the pair cannot be attributed to a simple Zeeman effect. 
Another curious feature is that the pair states were not 



observed'* in defect-free quantum dots with a physical di- 
ameter of 3000Aor equivalently a lateral confining length 
of 400A. 

The question then arises, what is the source of the ap- 
parent pair-tunneling peaks in the capacitance spectrum 
of the l//m quantum dots? As GaAs is not a supercon- 
ductor, it is surprising that a stable pair-tunneling state 
is observed at all. A theory of the elusive pair states 
in GaAs quantum dots must explain the following facts: 
1) the stability of a 2e~ tunneling state at zero magnetic 
field, 2) a pair instability at a critical field B = Bg, 3) the 
emergence of two distinct tunneling peaks for B > Be, 
and 4) the absence of pair-tunneling states in defect-free 
quantum dots with a confining diameter approximately 
< 400A. We propose here a model which is capable of 
explaining each of these observations. 

Although the magnetic field dependence might seem 
decoupled from the pairing mechanism, it turns out that 
an understanding of the field dependence is crucial to 
the formulation of a model for pair tunneling in quan- 
tum dots. Asymptotically, the energy of a parabolically- 
confined electronic state increases as huic/2 in a magnetic 
field. Here ujc is the cyclotron frequency. Inspection of 
the plots of Vg vs B in Fig. 6 of Ref. [lb] reveals that the 
energy of some of the tunneling states scales asymptoti- 
cally as huic/2, while others do not. States which exhibit 
the asymptotic huic/2 slope originate from 2D harmonic 
confinement on a length scale of 300A. Such confinement 
arises from the potential created by clustering of ionized 
Si-dopants in the AlGaAs backing layer*'®. Ashoori and 
co-workers have attributed the states whose energy is rel- 
atively insensitive to a magnetic field as arising from Si 
impurities that diffuse into the quantum well from the 
source lead*. In such hydrogenic impurities, the 1/r po- 
tential competes with the parabolic potential generated 
by the magnetic field to suppress the asymptotic huic/2 
slope. The effective Bohr radius of a Si impurity in GaAs 
is a*o K. lOOA. 

Consider now the pair states labeled X, Y , and Z . A 
consistent trend in these pair states is that above B^, 
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the two electronic states evolve with fundamentally dif- 
ferent slopes, indicating two distinct confining potentials. 
The upper state consistently has a slope close to 7}Wc/2 
whereas the slope of the lower state is indicative of that of 
a hydrogenic impurity. Consequently, the pair states ap- 
pear to be hybrid states composed of a Si-impurity and a 
parabolic well formed from the (5— dopants in the backing 
layer. A hybrid potential of this sort will give rise to two 
different length scales once the electron pair breaks up. 
Hence, such a potential will be the starting point for our 
analysis of the pair states. We point out that some exper- 
imental evidence exists in which both states of the pair 
evolve as 7}Wc/2 above Be, for example the data shown in 
Fig. 6a of Ref. [lb]. However, these states appear to be 
less stable and less reproducible than are the pair states 
of the X, Y, and Z type because they transform into the 
latter upon thermal cycling of the sample^. We take then 
the X, Y, and Z pairs to be the archetypal pair states 
and focus heretofore on explaining their origin. 

Barring something truly exotic, the only tool available 
to mediate pair tunneling is coupling to optical phonons. 
GaAs is a weakly-polar semiconductor with an optical 
phonon energy of 7}Wio=36.6 meV and coupling constant 
of a = 0.08. To determine whether phonons stabilize a 
pair state in a hybrid hydrogenic-parabolic potential on 
a quantum dot in an external magnetic field, we must 
solve the Hamiltonian H = -ffo(l) + Ho(2) + ?y/ri2 for 
the two and one-particle ground-state energies, E(2) and 
£"(!), respectively. In this expression®, rj = e^/eo with 
fo (=12.53) the static dielectric constant for GaAs and 
Ho(i) = Hi(i) + Hg-ph is the sum of a 1-body Hamilto- 
nian 

ff.C) = ^ - ^ + Vg(z) + if [(I - lo)' + r] 

+^(f)'(«^ + r) (1) 

and the electron-phonon interaction of the Frohlich type 
He-ph = y^ajafc 

k 

k 

where (P,r) are the electron momentum and position, 
g = (y^iray^^ , V is the crystal volume; the confining po- 
tential of the GaAs/AlGaAs quantum well of width L is 
Vq(z) = Vo 0(\z\ — L/2), ujo is the frequency of the oscilla- 
tion in the parabolic well formed by the (5-dopants and a|, 
creates an optical phonon with momentum k. In writing 
Eqs. (1-2) we have set k = m* = uj^o = 1- The coordi- 
nate ro = (xo,0,0) determines the distance between the 
Coulomb and parabolic potentials. We have also chosen 
the symmetric gauge A = ^(—y,x,0); considering the 
range of magnetic fields involved and based on an anal- 
ysis of D~ centers^'® in GaAs, we ignore paramagnetic 



contributions to Hi. Phonons will mediate a negative-U 
center if E(2) — 2E(l) < 0. Of course, this pair con- 
straint necessitates solving the coupled electron-phonon 
problem exactly. The relevant quantity that must be 
calculated is the bi-polaron attractive interaction, Eb- 
However, as GaAs is in the weak electron-phonon cou- 
pling regime, a = .08 << 1, it can be shown^ that the 
result based on perturbation theory-'^'^, Eb = 2Ep, is an 
excellent approximation, where Ep is the polaron bind- 
ing energy. Our approach will be to solve the electron 
problem separately using Quantum Monte Carlo (QMC) 
techniques. To this end, we redefine our working Hamil- 
tonian to be if = Hi(l) + Hi(2) + ry/ri2. Let E(2) and 
£"(1) be the exact ground state energies of iJ and Hi, re- 
spectively. The pair-binding condition can now be recast 
as E{2) + 2Ep - 2E{1) < 0. 

To obtain Ep, we need to solve for the phonon- 
enhancement in the binding energy of the ground state of 
Ho = Hi + Hg-ph- For the 3D free electron problem, the 
polaron energy is simply the Frohlich^^ result Ep = —a 
(in units of Tibj^o)- Platzman^^ has shown that for an 
electron bound to a Coulomb potential, the polaron bind- 
ing is enhanced in the strong coupling regime relative to 
the free polaron problem. We will show here that the 
polaron energy is further enhanced in our hybrid poten- 
tial even in the weak coupling limit. This enhancement 
assists the formation of a pair state. To proceed, we use 
Feynman's approach^^ and obtain the effective Euclidean 
action in terms of the electron path, Vt, 




by integrating out the phonon degrees of freedom. The 
presence of a retarded electron-phonon interaction, 1/r, 
and quantum well potentials precludes an exact analytic 
solution to this problem. We are guaranteed, however, if 
we define a trial action St 



-f Jdtr^-fJ dt (4) 

that the total energy will satisfy the variational up- 
per bound, E = Et — lim^^oo P~^{Seff — St) > 
Eq. We have defined to be the sum of the first 
and last two terms in the effective action and Et = 
— lim^^oo l3~^ In f Vrt exp[ST]- The parameters, C, W, 
K, and K^, are to be determined variationally. The po- 
laron binding energy is given by Ep = E — E(a = C = 0). 
Computation of Ep is compounded by the lack of central 
symmetry of the confining potential. Only for the case 
in which xq = can the integrals be performed in closed 
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form. It turns out that for a bound polaron, the ex- 
plicit xq dependence for the separations considered here 
amounts only to a few percent correction to the polaron 
binding energy^. Hence, we will use only Ep(xo = 0) so- 
lutions. We evaluated the path integrals using Laplace 
transform techniques to solve the saddle point equations. 
As the details of such calculations are rather lengthy, we 
only present the results for Ep. Results for the mag- 
netic field dependence of the polaron binding energy are 
shown in Fig. 1. We find that Ep has a weak field de- 
pendence. At 10 T, the increase in \Ep\ is 0.05 meV, 
which is five times smaller than the Zeeman energy (0.25 
meV). Hence, Ep is relatively insensitive to the magnetic 
field. At 5 = 0, the confining potential has enhanced the 
polaron binding energy from Ep = —a to Ep = —1.1a. 
For GaAs, this corresponds to an enhancement of 0.3 
meV or equivalently an enhancement of 0.6 meV for the 
bi-polaron binding energy. We will see that this enhance- 
ment plays a role in the pair-formation condition. 

We now need to evaluate accurately the energies in- 
volving the electronic degrees of freedom. For the range 
of magnetic fields we are dealing with, the ground state 
of the 2-fermion system is a spin-singlet, which means 
that the configurational part of its wavefunction must 
be symmetric. We can solve exactly this boson prob- 
lem by using projector QMC methods which are based 
on the property that in the asymptotic imaginary time 
(r) limit, the Euclidean evolution operator acting on 
a parent state $t projects out the ground state $o: 
$0 cx lim^^oo exp[— r(_ff — Eg)] ^t, where Eg is a 
suitable trial energy which shifts the zero of the energy 
spectrum. We transform the time-dependent Schrodinger 
equation for $ in Euclidean time r to a master equa- 
tion for the importance-sampled distribution P(TZ, r) = 

$T(7^)$(7^,r)/||V$7^||^ 



dP(TZ,T) 
dr 



ViP(TZ,T) - Fi(TZ)P(TZ,T) 



(EL(TZ)-Eg)P(TZ,T) 



(5) 



and use stochastic random-walks in configuration space 
(TZ = (ri,r2)) to solve this equation. Fi(TZ) = Vj- ln$T 
is the dnft velocity and El(TZ) = is the local 

energy. The trial function 



$^(7?.) = exp 



ri2 



2 -I- ari2 

2 

]^(?i(zi)exp 



[gi(l)g2{2) + gi(2)g2(l)] 



1 + a^r-i 



liPi 



(6) 



is used to guide the random- walk. In Eq. 6 gi(i) = 
exp[-&r?/(l + en)], g2{i) = exp[-iy{{x - xof + y'-')], 
(l)(z) = cos(kz) or cos(k_L/2) exp[A(_L/2 — \z\)] inside or 
outside the well, respectively, and with a, h, c, ji, v, ai, k, X 
variational parameters^. At sufficiently long times 
P(TZ,T co) $t(7?.) $o(^) (up to a normalization 



constant), where $o is the exact (nodeless) lowest energy 
state. In order to get this stationary distribution. Eg 
must be adjusted to be equal to the ground state energy 
E(2), given in turn by E(2) = lim^^oo {El{TI)) p(n,T)- 
As long as $t satisfies the right permutation symmetry, 
its functional form affects only the convergence and sta- 
tistical fiuctuations of £"(2) . In a similar way we calculate 
E{1). 

We now combine our variational estimates for the po- 
laron energy with our exact QMC calculations for £"(2) 
and £"(1) to solve the pair-binding condition, Ap = 
£(2) + 2Ep - 2£(1) < 0. When pair-binding occurs, 
Ap < and both electrons reside in the same electronic 
state with energy e = (£(2) -|- 4£p)/2 as measured in 
SECS. If Ap > 0, two distinct states are occupied: a 
lower one with energy ei = £(1) -|- £p and a higher 
one with energy £2 = E(2) + 3Ep — £(1). Using the 
method outlined above, we have calculated Ap at three 
separations between the hydrogenic and parabolic wells: 
xq = 5ao, xq = Soq, and xq = IOoq. We have plotted in 
Fig. 2 the energy of two electrons on the quantum dot 
as a function of the magnetic field. The salient features 
of these graphs are as follows: 1) at 5 = 0, pair states 
form only when xq exceeds a critical value Sag, 2) a 
magnetic field inhibits pair-binding at critical fields of 
0.2T (xq = Sflg) and 1.5T (xq = lOog) and 3) as in the 
experiments"'^, two new states emerge above Be; ei has 
a weak field dependence and £2 scales asymptotically as 
huic/2. The absence of pair-binding for xq < Sag is a 
signature that the electronic wavefunctions are too local- 
ized for the electron-phonon attraction to outweigh the 
Coulomb repulsion. Similarly, increasing the magnetic 
field further confines the wavefunction of the electrons 
in the pair state. As a result, the Coulomb repulsions 
increase and at a critical value of the field, the polaron 
effect can no longer hold the pair together, resulting in 
pair-breaking. 

Of course, the confining potential caused by the 6- 
dopants is only approximately parabolic. Anharmonic- 
ities are undoubtedly present. We explored the role of 
an anharmonic correction of the form ( ((x — xq)^ + y'^)^ . 
This correction increases the barrier between the hydro- 
genic and parabolic potentials (see inset in Fig. 1) thus 
increasing the mean distance between the two electrons 
and consequently reducing the Coulomb repulsion. The 
results shown in Fig. 3 are for a separation of xq = Sag. 
As is evident, anharmonicity of this type increases pair- 
binding and pushes Be to a higher value of 2T. In 
addition, the magnetic field dependence of ei and £2 is 
as in the harmonic case. We conclude then that the ex- 
perimental observations of pair-binding can be explained 
by a model in which optical phonons bind two electrons 
in a Si impurity and a roughly parabolic well (created 
by the (5-dopants) that are sufficiently far apart^'*. The 
critical separation even when anharmonic corrections are 
included seems to be Tag = 700A. We would predict 
then an absence of pair-binding in quantum dots that 
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are either clean (no Si) or laterally smaller than Tag in 
confining diameter, consistent with experiments reported 
in Ref. [lb]. We close by pointing out that the polaron 
binding considered here increases monotonically as the 
GaAs quantum well width diminishes, and in the strict 
2-dimensional limit Ep = — 7ra/2. Such an enhance- 
ment is sufficient to cause pair-binding in samples where 
both lateral confining radii and quantum well widths are 
reduced. Hence, a critical test of the model reported here 
would be a detailed sample-size dependent study of the 
stability of pair states in GaAs quantum dots. 
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FIG. 2. Electron energies as a function of magnetic field 
B. For each pair-potential separation xo, the lower branch 
is ci — Ep while the upper one corresponds to ^2 — Ep. The 
open circles represent the QMC results whose statistical er- 
ror bars are 0.06 meV (smaller than the size of the circles). 
The curves are spline interpolations to the QMC data. The 
dashed-dotted line represents the field dependence of the low- 
est Landau level fta;c/2. 

FIG. 3. Effect of anharmonicities on the electron energies 
at pair-potential separation xo = Sag. Note anharmonic cor- 
rections shift Be to higher values. The meaning of the curves 
and open circles is the same as in Fig. 2. 



4 



